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Bk from quenched QCD with exact chiral symmetry 
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We present a calculation of the standard model AS = 2 matrix element relevant to indirect CP 
violation in K — > nn decays which uses Neuberger's chiral formulation of lattice fermions. The 
computation is performed in the quenched approximation on a 16 3 x 32 lattice that has a lattice 
spacing a ~ 0.1 fm. The resulting bare matrix element is renormalized non-perturbatively. Our 
mam result is B% GI = 0.87(8)^1^, where the first error is statistical, the second is systematic and 
the third is an estimate of the uncertainty associated with the quenched approximation and with 
the fact that our kaons are composed of degenerate s and d quarks with masses ~ m s /2. 
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Introduction. — Ginsparg- Wilson (GW) fermions [1-4], 
with their exact SU(Nf)i l x SU(Nf)ji chiral-flavor sym- 
metry at finite lattice spacing [5], provide unique oppor- 
tunities for exploring the physics of light quarks through 
numerical simulations of lattice QCD (e.g. reviews in 
[6, 7]). In this regularization the AS = 1 effective weak 
Hamiltonian renormalizes with the same pattern as in 
the continuum. As a consequence, in the presence of an 
active charm quark, the AI =1/2 rule in K — » 7T7t de- 
cays can be studied from the simpler K — > it amplitudes 
without having to perform difficult power subtractions 
[8]. In addition, this exact chiral-flavor symmetry im- 
plies full O (a) improvement. In this letter we present 
the results of a quenched calculation of the matrix ele- 
ment of the AS — 2 effective weak Hamiltonian relevant 
to K°—K° mixing, performed using Neuberger's imple- 
mentation of GW fermions [3, 4]. The bare matrix ele- 
ment is renormalized non-perturbatively in the RI/MOM 
scheme a la [9]. This calculation is the first study of 
a matrix element of a four-quark operator using a lat- 
tice formulation which has exact chiral-flavor symmetry 
at finite lattice spacing. It calls upon many of the in- 
gredients required for studying the AI = 1/2 rule or e' 
(three-point functions, non-perturbative renormalization 
of four-quark operators, etc.), but avoids the numerically 
challenging computation of eye diagrams. Furthermore, 
because this AS = 2 matrix element has been extensively 
studied with other fermion formulations (e.g. reviews in 
[10, 11]), it provides a good test of the reliability of our 
approach. 

In addition to using Neuberger fermions, we make a 
number of other improvements on the methods tradition- 
ally used to compute AS* = 2 matrix elements. In order 
to eliminate sizeable unphysical finite-volume contribu- 
tions from topology-induced fermion zero-modes, we use 
the fact that the relevant AS = 2 operator is purely left- 
handed and consider correlation functions constructed 
from left-left quark propagators only. These propaga- 
tors are free from zero-mode contributions, as the lat- 



ter are chiral. Thus, to create and destroy kaons at 
rest, we use the time component of left-handed, quark- 
bilinear currents. The use of left-handed currents was 
proposed in [12, 13] in the context of the e-regime of 
Gasser and Leutwyler [14]. We also use chiral sources in 
the RI/MOM non-perturbative renormalization (NPR) 
of our operators to avoid zero-mode contributions, which 
would otherwise appear. This has the added advantage 
of eliminating the leading chiral-symmetry-breaking con- 
tributions to the NPR Green functions. 

K°—K° mixing induces indirect CP violation in K — > 
7T7T decays, quantified by the parameter e for which the 
theoretical expression [15] and experimental value [16] 



ee"** ~ C e CK-(//)B*MA 4 A 10 7j[(l - p) m S (x t ) (1) 

+P{x t ,x c , ...)] = (2.282 ± 0.017) x 10" 3 , 
with B K defined through (F w = 93 MeV) 

1 fi 

(K \O As ^)\K ) = -M\Fl x B K 0*) (2) 



and 



O as =2 = [*y„(l - 75)dp7 M (l - 7 5 )d] • (3) 



Here, C e is obtained from well measured quantities, A, A, 
p and fj are Wolfenstein parameters [17] and 772, Ck, Sq 
and P incorporate perturbative, short-distance physics 
(P also contains CKM factors). Combined with a de- 
termination of Bk-i the measurement of e provides an 
important hyperbolic constraint on the summit (p, fj) of 
the unitarity triangle, as indicated by Eq. (1). It is ad- 
vantageous to compute Bk instead of the AS = 2 matrix 
element, because many lattice systematic errors cancel in 
the ratio that defines Bk- A preliminary version of the 
present calculation was reported in [18]. Preliminary re- 
sults obtained by the MILC collaboration with a different 
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implementation of GW fermions were also published in 
the same volume [19]. 

Computational details. — The simulation is performed 
in quenched QCD with (3 = 6.0 and V = 16 3 x 
32. We use a sample of 80 gauge configurations, gen- 
erated with the standard Wilson gluonic action, re- 
trieved from the repository at the "Gauge Connec- 
tion" (cf. http://qcd.nersc.gov). Fermion propaga- 
tors are obtained from a local source, using Ncu- 
berger's action [3, 4] with bare quark masses am = 
0.040,0.055,0.070,0.085,0.100. The sign function of the 
Hermitian Wilson-Dirac operator, X, is obtained by 
making an optimal rational approximation [20-22], after 
explicit evaluation of the contributions from the lowest 
eigenvectors of X^X. The computation of the propaga- 
tors thus uses nested multi-conjugate gradient inversions 
(for more details see [23, 24]). The lattice spacing is de- 
termined from r = 0.5 fm, where r is the Sommer scale 
[25]. This gives a" 1 = 2.12 GeV. 

Meson correlation functions and fits. — From the quark 
propagators we compute the two-point function 

Cjj(xo) =Y,(Jo( x )M0)) , (4) 

X 

and the three-point function 

Cj O j(x ,y ) = J2( J ^ x )°As=2(0)My)) , (5) 

where J M = ^(1 - 7 5 )s, q = (1 - jrD)q, p = 1.4 (see 



d and 0As=2 



[23, 24]), J is obtained from J with s 
is obtained from Eq. (2) with d — > d. Statistical errors 
are estimated with the jackknife method. 
The ratio 



R(x ,yo) = o 



3 Cj j(x ,yo) 



8Cjj(x )Cjj(y ) 



(6) 



is fitted to the asymptotic form (a <C x <C yo "C T): 

R(x ,y ) -^B h K ^ , (7) 



where the fit parameter B^? rc is the bare bag parameter. 

The kaon mass Mk, on the other hand, is determined 
from a fit of the time-symmetrized, two-point function to 
the standard asymptotic form (a <C x ) 

Cjj(x ) — > ^e- M - T ' 2 cosh [M K (T/2 - x Q )} , (8) 

with Z = \(K°\J o \0)\. We find that the two- 
point function is asymptotic from ax = 5 on and 
that the statistically optimal fitting range for aMx 
is 5 < axo < 12, 13, 14, 15 and 16 when am = 
0.040, 0.055, 0.070, 0.085 and 0.100, respectively. 

For the ratio of Eq. (6), we find that the asymptotic 
behavior sets in for oxq > 5 and ayo < 27. Because our 



lattice is rather short in the time direction, we have to 
worry about time-reversed contributions. Assuming that 
the time-reversed matrix element is approximately the 
same as the forward one and that the two-kaon energy 
is approximately 2Mk (i.e. that finite-volume effects are 
not significant), then for ax < Q,ay > 26 the time- 
reversed contributions never exceed 1.5% of the forward 
signal. We also checked this explicitly by a fit which 
includes time-reversed contributions and allows for finite- 
volume shifts on the time-reversed matrix element and 
the two-kaon energy. We thus use the range 5 < ax < 6 
and 26 < ay < 27 to calculate our observables. All 
fits are excellent and our results for aM K and B^ rc are 
summarized in Table I. 



aM K 



= (a) 



0.040 0.253(5) 
0.055 0.288(4) 
0.070 0.321(3) 
0.085 0.352(3) 
0.100 0.382(2) 



0.70(7) 
0.71(6) 
0.73(5) 
0.74(4) 
0.75(4) 



TABLE I: Mesons masses and B^ arc (a) vs quark mass. 



Non-perturbative renormalization. — We perform all 
rcnormalizations non-perturbatively in the RI/MOM 
scheme a la [9]. Thus, we fix gluon configurations to 
Landau gauge and numerically compute appropriate, am- 
putated forward quark Green functions with legs of mo- 
mentum p = \fp*. We define the ratio 



TZ RI (m,p',g ) 



Tj(m,p 2 ) 2 



(9) 



where To is the value of the non-perturbative, amputated 
Green function of operator O computed with left-handed 
quark sources and projected onto the spin-color structure 
of O. We compute this ratio for all five quark masses and 
perform a linear chiral extrapolation in m 2 as prescribed 
by a next-to-leading order expansion of the Green func- 
tions in quark mass. We find that the mass dependence 
is very mild and is well described by this linear form. We 
must then isolate the "perturbative part" of this ratio to 
get the renormalization constant appropriate for renor- 
malizing B^ rc (a). A straightforward operator product 
expansion (OPE) in 1 jp 2 yields the following p 2 behav- 
ior for lZ m : 



n RI (0,p z ,g o ) 



4 + ^;; J (<yo)xp |< V) 



*>RGI 



RI/ 



-B x (ap) 2 



(10) 



where Z R( ^ l (go) is the renormalization constant that 
takes B^ arc (a) to the renormalization-group invariant pa- 
rameter -B]^ GI , and U Rl (p 2 ) describes the running in 
p of the corresponding Rl-schcmc renormalization con- 
stant Z Rl (p 2 ,go)- The ellipses in Eq. (10) correspond 
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to higher-order terms in the OPE and higher-order dis- 
cretization errors. 

To describe the running of Zg l K (p 2 , go), we use the 2- 
loop expression obtained by combining the MS— NDR re- 
sult of [26] and the MS-NDR -> RI matching result of 
[27], with N f = 0. Thus, 



U™(p*) = a s (p 2 ) 



2 < U 



Pi a s (p 2 ) 



t? 1 



70 



(11) 



with /3 = 11, /3i = 102, 70 = 4 and jf 1 = 287/3 - 
176 In 2. The strong coupling constant that we use is 
taken from [28] and corresponds to a s (fj^ ef ) — 0.0926(12) 
in the MS scheme with [i m { — 94.1(3.6)/ro. To obtain the 
coupling at other scales, we integrate the 2-loop running 
equation exactly and solve it numerically. 

Our results for lZ m (0 , p 2 , go) are plotted in Fig. 1 as 



a function of p 2 



Eu=o( sina 2V) 2 - Usin S the 



lattice definition of momentum significantly reduces dis- 
cretization errors. Also shown is the fit of these re- 
sults to the functional form given by Eqs. (10) and (11) 
with Zg^ l (go), A and B as parameters, in the range of 
2GeV 2 < p 2 < 10GeV 2 . We find that the OPE and 
discretization error terms kept in Eq. (10) are sufficient 
to describe the data in this range. The fit actually de- 
scribes the data in a much larger range, indicating that 
the retained terms dominate. The results of the fits are 
summarized in Table II, where we also display a fit to the 
extended momentum range 0.65 GeV 2 <p 2 < 15.9 GeV 2 
with additional 1/p 4 and (ap) 4 terms, and a fit us- 
ing the continuum p 2 in the range 2 GeV 2 < p 2 < 
10 GeV 2 . All fits are excellent and produce compatible 
results. The renormalization constants at 4 GeV 2 in the 
RI and MS— NDR schemes are obtained by multiplying 
Zjj^'(jo) by the appropriate two- loop running expres- 
sions £/ RI (4GeV 2 ) and C/ NDR (4 GeV 2 ), respectively, with 
Nf and a s chosen as above. £/ NDR (p 2 ) is given by Eq. 
(11) with 7 RI -> 7^ DR = -7 for Nf = 0. 

Physical results. — The central values for B K iIO (a), 
Zlf{9o), ZNDR (4G eV 2 , 3o ) and Zf K (4GeV 2 , 5o ) are 
obtained as described above. While statistical and sys- 
tematic errors in the B-parameters and renormalization 
constants will be correlated in our final results, we also 
wish to give results for the renormalization constants 
themselves. The systematic errors on these constants 
take into account the following sources of uncertainty: 
the errors on a s (p 2 cl ) and /x rc f, given above; the varia- 
tion due to a ±10% uncertainty on the lattice spacing, 
which is typical in quenched calculations [29]; the dif- 
ference between the renormalization constants obtained 
using the lattice and continuum definitions of quark mo- 
mentum. The latter yields an estimate of discretization 
errors which turns out to be the dominant uncertainty. 
We take it to be a symmetric error. Adding all of these 
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FIG. 1: Tl m (0,p 2 ,g ) vsp 2 . The solid triangles correspond 
to the points used in the fit to Eq. (10) while the open ones 
are not included. The curve is the result of the fit. 



variations in quadrature, we obtain: 



(4 GeV 2 , g ) 
Z R ^(4GeV 2 , 50 ) 



^NDR 



1.261(9)^ 
0.908(6)^ 
0.897(6)1? , 



(12) 



for p 



1.4. For comparison, we compute the 



renormalization constant Z R ^ PT at one-loop in per- 
turbation theory, matching the results of [30] with 
the standard RI/MOM scheme used here. We find 
Z R |; PT (4GeV 2 ,3o) = 0.95(5) with the value of a s given 
above. The central value is obtained by expanding the 
ratio of ^^2(4 GeV 2 ,. g ) to Z PT ( ff0 ) 2 in a s . The un- 
certainty is determined from the errors on a s (/i 2 of ) and 
/i rG f and from an estimate of O (a 2 ) corrections obtained 
by considering the unexpanded ratio. While the one- 
loop correction to Z A g =2 (4GeV ,go) is very large (of 
order 70%) and the perturbative result for ^ PT (go) 2 
is approximately 15% below the non-perturbative value 
found in [23], the agreement of Z RI,PT (4 GeV 2 , go) with 
the non-perturbative value given in Eq. (12) is surpris- 
ingly good. This is due to the large cancellations in 
^A5=I(4GeV 2 , 5o )/^ T (3o) 2 . 

We now turn to Bk- Our lightest pseudoscalar meson 
is very close to having the mass of the kaon. Since our 
results for Bk are linear in M K , we extrapolate them lin- 
early to the physical point, as shown in Fig. 2. To extrap- 
olate to lighter quarks, chiral logarithms should be taken 
into account. Thus, we also perform a fit to the func- 
tional dependence predicted by one-loop quenched chi- 
ral perturbation theory (Q^PT) [31], supplemented by 
an (M/AttF) 4 term to parametrize higher-order contri- 
butions, which are expected for our pseudoscalar meson 
masses. Here, F is the leading-order leptonic decay con- 
stant, which we take to be F n , and M is the leading-order 
pseudoscalar meson mass, which we set equal to meson 
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p^-range [GeV 2 ] (p 2 def.) 




W [GeV^] 


{apY 


1/p 4 [GeV 4 ] 


(apf 


[2, 10] (lattice) 
[0.65, 15.9] (lattice) 
[2, 10] (continuum) 


1.261(9) 
1.264(13) 
1.269(8) 


-0.042(16) 
-0.035(23) 
-0.051(16) 


0.028(2) 
0.024(5) 
0.022(2) 




-0.013(13) 





0.0010(9) 




TABLE II: Rernomalization constants and OPE coefficients obtained from fits of 7?, RI (0, p 2 ,go) to the OPE form of Eq. (10) 
for different momentum ranges and for the lattice and continuum definitions of p 2 . 



masses obtained in our simulation. The fit parameters 
are the RGI value of the chiral- limit i?-parameter, _B RGI , 
the scale of the one-loop logarithm and the coefficient of 
the (M/47rF) 4 term. The fit, also shown in Fig. 2, is 
only meant to indicate how a chiral extrapolation could 
lead to much smaller values of the £>-parameter in the 
chiral limit, such as the one found in the large- N c cal- 
culation of [32]. We would need results at lower masses 
and with better statistics to confirm the presence of the 
logarithm. For the time being we find, in the quenched 
approximation, B RGI = 0.53(11) j^Q , where the central 
value is obtained from the xPT fit- The first error here 
is statistical, the second results from implementing the 
systematic variations that were discussed above and the 
third is the difference between the central values of linear 
and xPT fits. 

As evident from Fig. 2, the value of B K Gl at the phys- 
ical point is insensitive to the choice of functional form 
in the fit. We thus use the simpler linear fit to determine 
B RGl at that point. We allow for a 20% uncertainty in 
the determination of the strange quark mass, which is 
typical of the variations observed in quenched calcula- 
tions of this quantity [29]. In addition, we account in a 
correlated fashion for the statistical errors and the sys- 
tematic variations that were discussed in the paragraph 
on the renormalization constants. We do not include a 
discretization error on the bare value of Bk, because we 
have no way of estimating it from our simulation per- 
formed at a single value of the lattice spacing. However, 
-B-parameters are ratios of very similar matrix elements, 
in which some discretization errors should cancel. More- 
over, the limited experience that we have with Neuberger 
fermions suggests that discretization errors on quantities 
such as Fk are small at (3 = 6.0 [33]. 

Our results for Bk in various schemes are summarized 
in Table III. They are in excellent agreement with world 
averages for this quantity (e.g. [10, 11]), which are based 
on quenched, staggered results [34-36]. However, our 
precision is not sufficient to exclude the rather low values 
found with domain- wall fermions [37, 38] . 

In addition to the systematic errors accounted for in 
the results of Table III, one must also consider the error 
associated with the fact that our kaon is composed of 
degenerate quarks with masses ~ m s /2, instead of an s 
and a d quark. This error is thought to be roughly 5% 
on the basis of xPT estimates reviewed in [39]. In this 
same review, Sharpe suggests adding a 15% quenching 




M| [GeV 2 ] 



FIG. 2: The renormalization group invariant parameter 
as a function of kaon mass squared. The solid circles are the 
result of our simulation which are fitted to a line (solid curve) 
and the QxPT expression (dashed curve). The solid triangle 
is the value of -B RGI at the physical point. The squares cor- 
respond to the chirally extrapolated values _B RGI as obtained 
from the linear (solid square) and QxPT (open square) fits. 



RGI NDR (4 GeV 2 ) RI (4 GeV 7 )' 
B K 0.87(8)+i 0-63(6)+j 0.62(6)+^ 



TABLE III: Final results for Bk at the physical point and in 
the chiral limit in different schemes. The first error is statis- 
tical and the second is systematic (see text for a discussion of 
additional m s = rrid and quenching systematics) . 



uncertainty. We have accounted for part of this error in 
considering the variations due to the uncertainties in the 
determination of the lattice spacing and of the strange 
quark mass. We find these variations to be very much 
smaller than 15%, which might suggest that the quench- 
ing error estimate given in [39] is rather conservative. 
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